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We develop a unified theory for clocks and gravimeters using the interferences of multiple atomic 
waves put in levitation by traveling light pulses. Inspired by optical methods, we exhibit a prop- 
agation invariant, which enables to derive analytically the wave function of the sample scattering 
on the light pulse sequence. A complete characterization of the device sensitivity with respect to 
frequency or to acceleration measurements is obtained. These results agree with previous numerical 
simulations and confirm the conjecture of sensitivity improvement through multiple atomic wave 
interferences. A realistic experimental implementation for such clock architecture is discussed. 



Atom interferometers p] have been used successfully to 
build accurate optical clocks [5] and gravimeters [21 H] ■ 
However, in the presence of a gravity field, these systems 
faced critical limitations caused by the free fall of the 
atoms used in the interferometric measurement. Among 
these, we mention the finite time of atomic interroga- 
tion, and the curvature shifts resulting from the trans- 
verse field exploration by the traveling atoms [5]. 

Recently, we have exposed a strategy to overcome 
these limits, which relies on the atomic levitation thanks 
to a series of vertical traveling light pulses [HI [3- ^ 
unique laser field performs simultaneously the atomic 
confinement and the atomic interrogation, which are 
usually operated by distinct laser fields in typical 
setups |8j. In our initial proposal [6 , the setup jug- 
gles with a single atomic cloud thanks to a series of 
synchronized velocity-sensitive Raman 7r-pulses. The 
feasibility of this strategy was demonstrated by Hughes 
et al. [5j in a similar experiment achieving over a hun- 
dred bounces. In a second proposal [7], we considered 
a different atomic sustentation process, which exploits 
instead the quantum interferences of multiple flying 
wave-packets. These packets bounce on a series of short 
7r/2-pulses, which give rise to a myriad of trajectories 
and yield a multiple-wave atom interferometer |10j in 
levitation. The phases of the light pulses are finely 
tuned in order to induce constructive interferences in 
a network of preferred levitating paths. An analog of 
such quantum trampoline has been realized recently by 
Robert-de-Saint- Vincent ct al. [TT| . 



preparing the way towards the implementation of such 
clock architectures. 

We proceed as follows. Sec|l] presents the theory of 
the atomic phase-matching with two separated tt/2 light 
pulses. We expose in greater detail the argument of 
Ref. if and derive the propagation invariant used in the 
subsequent developments. Sec|ll] recalls the principle of 
the atomic sensor using multiple- wave atomic levitation. 



In Sec III we derive the levitating wave-function when 



the system operates at resonance. In Sec IV we establish 
the levitating wave-function when the system operates 
out of resonance, and obtain the sensor sensitivity as a 
clock or as a gravimeter in the short-pulse regime. We ob- 
tain analytical expressions for the fundamental sensitiv- 
ity limit of such systems, in remarkable agreement with 
previous numerical simulations presented in [HI [7]. In 
Sec|Vj we exhibit atomic systems suitable for this setup 
with realistic experimental parameters. 



I. ATOMIC PHASE-MATCHING WITH A 
SEQUENCE OF 7r/2-PULSES. 

We consider a dilute two-level atomic sample well- 
described by a two-component wave-function, noted 
{'ipb{z,t),il)a{z,t)), and initially entirely in one of these 
two energy levels. This sample experiences a free fall in 
the uniform gravity field, thus following the Hamiltonian 



A moving two-level atom can be considered as a clock, 
and significant advances can be achieved by adopting a 
unified point of view for clock and inertial sensors |12j . 
This is exactly the philosophy of our approach, which 
gives a complete quantum description of the levitating 
atoms. Preliminary steps towards the construction of 
a levitating atomic gravimeter have been taken [9l lllj . 
In contrast, levitating atomic clocks still await their 
experimental realization. The purpose of this article 
is to explore in greater detail the behavior and the 
sensitivity of such suspended atomic systems. We 
also bring out suitable atomic species for these setups. 
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A key idea of the levitating multiple-wave interferometer 
is the fine-tuning of atomic phases through a sequence 
of 7r/2-pulse pairs. In order to expose this concept, it 
is useful to give a short reminder on the propagation of 
ideal atomic waves in time-dependent quadratic poten- 
tials with the ABCD method [H]. Since the external po- 
tential is one-dimensional, we restrict without loss of gen- 
erality the discussion to ID wave-packets, taking an ini- 
tial sample wave-function %lji{z,tQ) = wpi(z, Zca,Pca,wa) 
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of the form 



Wpi{z,ZcO,PcO,Wo) = 



(2) 



Using the ABCD theorem for atomic waves [12], one can 
express this wave-packet at a later instant t > as 

Mz,t) = wpi{z,zS),Pc{t),w{t))e'''^^'^'"P^<>-''-'«^/^ (3) 

The central position Zc{t) and central momentum Pc{t) of 
the wave-packet follow the classical equations of motion 
of a point-like particle evolving in the considered Hamil- 
tonian, situated at the position Zco with the momentum 
Pco at the instant t = Iq. The term Si{zcQ,Pco,t — to) is 
the action of the particle in the internal state integrated 
along the classical trajectory 



Si{z,o,PcO,t-to)^ f dt' - Ei-mgz,{t') 

Jto 



2m 



The width w is simply given by 



(4) 
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It is useful to express the ABCD theorem in the momen- 
tum picture. The Fourier transform of the wave-packet 
of Eq. ([2]) gives readily 



V'i(p,*o) 
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Performing a Fourier transform on the packet of Eq. ^ , 
one immediately obtains the momentum wave-packet at 
a later instant 



(p-Pcr 

2mg ^lSi(ZaO,PcO,t-to)/h^~lZap/h 



. . (7) 

We have introduced the initial and final width in mo- 
mentum Wpo = h/wo and Wp = h/w. The effect of a 
short 7r/2-pulse of wave vector k = efcz [i] on an atomic 
wave-packet is efficiently described by a 2 x 2 Rabi matrix 
evaluated at the central instant tc of the pulse 



1 / I j g-i(efcz-a)t) 
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In this expression £ is a position operator. The non- 
diagonal matrix elements are thus quantum operators 
involving the laser phase (l){z,t) = ekz — cot, with e — I 



[i] We adopt the convention that fc > is the absolute value of the 
wave- vector associated with all the considered pulses, and we use 
e = 1 and e = — 1 for an upward and a downward-traveling pulse 
respectively - the axis Oz is oriented upwards -. 



and e = — 1 for an upward and downward-traveling pulse 
respectively. Their action onto an atomic wave-packet 
reflects the acquisition of a laser phase 0(zc,ic) a-nd of a 
momentum quantum through the atomic recoil 
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X wpi{z, ZcPc + £fik,w) . 



The proposed experiment rests on the ability to perform 
a controlled momentum transfer from the light field to 
the atoms thanks to a sequence of 7r/2-pulses. In this 
view, it is important to understand how such transfer 
can be done with two 7r/2-pulses separated by a time 
interval. There is an obvious similarity with the Ramsey 
interrogation, but the novelty of our approach is that 
it takes into account the quantized atomic motion in the 
gravitational field. Fig{T]illustrates the effect of this pulse 
sequence on an atomic sample. We consider two short 




Time 



FIG. 1: Action of a pair of copropagating 7r/2-pulses on a free- 
falling atomic wave-packet. The phase difference between the 
two outgoing wave-packets comes from the classical action 
and laser phase acquired on each path, and from the distance 
of their centers. 

upward 7r/2-pulses performed respectively at the times 
U = and tf = T on the sample and with the phase 
expressed with the central time as reference: 

= kz-u:i{t~T/2), (j)2iz,t) ^ kz-Ljj2{t-T /2)+(t)c 

(10) 

Note that the evaluation of the phase at the central time 
is an exact result, which follows from the ttt theorem |13j . 
and not a mere approximation. The gravitational accel- 
eration yields different resonance conditions for the first 
and the second light pulse, to which the frequencies wi 
and UJ2 are adjusted respectively. The pulses split the 
atomic cloud into four wave-packets, following two in- 
termediate paths labeled with the internal atomic states 
a, b. The upper-state and lower-state wave-functions re- 
ceive a contribution from each path, corresponding to 
wave-packets of different central altitude but of common 
central momentum. In order to describe the quantum 
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interferences between the atomic wave-packets after this 
sequence, it is thus easier to work in the momentum pic- 
ture, for which the packets of a given energy share the 
same center. We consider an initial lower-state wave- 
packet 








(11) 



of a similar form as Eq.Q up to a different constant 
phase. This choice will turn out convenient to compute 
the levitating wave-function. Using the Rabi matrix of 
Eq.Q, and the ABCD law in momentum picture given 
by the Eqs.([6|7|, one readily obtains the form of the 
wave-function after two 7r/2-pulses 



(P-Pcb) 



El. 
El. 



(12) 

We have noted Zfa, Pfa = Po ^ mgT and Zfh,pf(, — 
Po + hk — mgT) the respective central position and mo- 
mentum of the lower-state and upper-state atomic packet 
immediately before the second light pulse [See Fig T] . In 
the equation above, the phase $ei.e2 is associated with 
the wave-packet that was in the internal state ei after the 
first light pulse and in the internal state 62 after the sec- 
ond pulse. We note 0ei,e2 the contribution to this phase 
arising from the atom-light interaction. The expression 
of the phases $ei,e2 follows from the ABCD theorem in 
momentum picture and from the Rabi matrix 



(13) 



We have used the short-hand notations Sa = 
Sa{zo,Po,T) and Sb = Sbizo.po + hk,T). The phase 

(^61,62 read 0a, a = 0, (t)aM = 4>2{zfa,T), (t)b,a = (l>l{zo,Q)- 

(f>2{zfb,T) and 4>b,b = 0i(2^o,O). In order to obtain a full 
momentum transfer for the sample, one aims at produc- 
ing constructive interferences in the upper atomic state 
and destructive interferences in the lower atomic state, 
i.e. the outgoing wave-function ^p{p, T) should be of the 
form 



^ 



/ 

(14) 

This is possible only if the phases ( 13 1 satisfy the rela- 
tions 



V 



Zfb{p-po) 



*a,6 = * 
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(15) 



In order to preserve the norm of the double component 
wave-function ip during the unitary evolution, these rela- 
tions are not independent: it is for instance sufficient to 
ensure constructive interferences in the upper state in or- 
der to obtain destructive interferences in the lower state. 



We thus choose without loss of generality to satisfy the 
condition ^a,b = $6,6. Considering the phase reference 
used for the pulses [see Eq.(lO)], one sees that the pre- 
vious relation is satisfied if and only (j)c — and if the 
pulse frequencies are adjusted to the values 



<2 
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with pi = pq and P2 = Po ^ mgT. These equations state 
simply the energy conservation for each pulse, taking into 
account both the change of internal level and the atomic 
recoil. One should note that the phase matching given by 
the relations (16 1 is not sufficient to guarantee efficient 
constructive interferences in the excited state. One must 
also require a sufficient overlap of the interfering wave- 
packets, namely Wp{zfb — Zfa) 1, or equivalently 



^'^Dopplcr^ « 1 



(17) 



with Awj^pppjg^. — kwp/m the Doppler frequency spread 

experienced by the atomic sample because of its dis- 
persion in momentum. We note that this condition 
can also be expressed in terms of Ramsey interrogation 
as AwDoppier < Acjj^amsey/l^'r) since Awj^^msey = 
2tt/T is the frequency resolution expected from the Ram- 
sey interrogation with the considered sequence of tt/2- 
pulses. If this condition is not met, the momentum 
spread reduces the constructive interferences through a 
partial overlap of the atomic wave-packets. Alternatively, 
this condition can be stated as w 3> VrT: the vertical co- 
herence length of the atomic wave-packets must be much 
greater than the separation of their respective centers in 
order to yield constructive interferences. This is thus a 
perfectly intuitive result. 

We establish now a key property: the phase ^a,b = ^b.b 
picked up by the wave-packet ending in the upper state 
is not only path-independent but also independent from 
the initial position zq, providing us with an invariant [i] 

I = PiZi — pfZf + S + hcj) . (18) 

To show this, we compute the quantity / on each path, 

Iaizi,pi,T) = piZ^-pf{zi+piT/m- gT'^/2) (19) 

+ Sa{zi,Pt,T) + 02(^4 +piT/m - gT^,T) 

Ib{zi,Pt,T) = piZi~pf{zi + {pi + hk)T/m~ gT^/2) 

+ Sb{z,,p, + hk,T) + (j}^{z,,Qi), (20) 



when the resonance conditions of Eq. ( 16 1 and the condi- 
tion 0c = are satisfied. A straightforward computation 



[i] We correc t he re a typal in [7] , which explains the difference be- 
tween Eq.(|18[l and the invariant considered in [7]. 
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gives Ia{zi,pi,T) = Ib{zi,pi,T) = I{pi,k,T) with 



I{k,p„T) 



2 1 



(Eb + Ea) 



P 

2m 
1 



T 



hk)gT^ - -mg^T^ (21) 



At a deeper level, the path-invariance of the quantity 
/ is a consequence of an important feature of atomic- 
wave propagation, which becomes obvious in a 5D- 
formalism [T3]: the center of the wave-packet follows an 
equiphase trajectory with a generalized phase including 
an extra-degree of freedom, namely the proper time con- 
jugated to the atomic mass. As announced, the invari- 
ant / depends on the initial momentum but not on the 
initial position. This property allows one to address si- 
multaneously the levitating atomic packets situated at 
different altitudes with a single laser frequency. Besides, 
this fact turns the momentum picture into a very con- 
venient framework for the description of the levitating 
interferometer. 



II. PRINCIPLE OF THE MULTIPLE- WAVE 
ATOMIC LEVITATION. 



Altitude 




FIG. 2: Levitating atomic trajectories in the sequence of 
pulses. The first four pulses generate a vertical Borde-Ramsey 
interferometer. The central positions of the wave-packets ex- 
plore a network of paths whose number doubles at each laser 
pulse. 



For convenience, we expose in this Section the princi- 
ple of the sensor using multiple- wave levitation, following 
the lines of Ref. ^7 . In order to obtain the desired sus- 
pension of the atomic sample, one extends the previous 
pulse sequence as follows. Immediately after the pair 
of upward-travelling 7r/2-pulses, one performs a pair of 
downward-travelling 7r/2-pulses with an identical inter- 
pulse duration. The four-pulse sequence obtained consti- 
tutes a Borde-Ramsey interferometer, bent here by the 
gravitational field. This interferometer corresponds to 
the trajectories in the dashed box of Fig. [2j The previ- 
ous discussion shows that this second sequence performs 
a nearly full population transfer if the energy conserva- 
tion is verified for each pulse - using again the mid-time 
of the pulse pair as phase reference - and if the momen- 



tum dispersion is small enough as to fulfill Eq.(17|. One 



then achieves a net momentum transfer of Ap = 2hk per 
atom. Besides, for the special interferometer duration 



hk 

mg ' 



(22) 



the atoms end up with their initial momentum value. 
For this special duration, by repeating the four-pulse se- 
quence, one obtains a periodic atomic motion in momen- 
tum, thereby enabling levitation if the sample is initially 
at rest. Thanks to this periodicity, the energy conser- 
vation relations ( 16 ) are identical for each four-pulse se- 
quence. Indeed, these conditions yield only two differ- 
ent resonant frequencies splitted from the transition fre- 
quency UJab = [Eb — Ea)/h by the recoil frequency 



he 

2m ' 



OJ- 



2,3 



he 

2m 



(23) 



The first and the fourth pulse, as well as the second and 
the third ones, have identical resonant frequencies. 

When the previous resonance conditions (T :— 
2o,t^i.2,3,4 := "^1234) ^'^'2 fulfilled, and for a sufficiently 



narrow atomic momentum distribution (17), the atomic 



sample spreads through the network of suspended trajec- 
tories sketched on Fig. [2] As explained in Ref. [7] , when 
the parameters (T, a;i^2.3,4) differ from these resonant val- 
ues, the phase-matching condition is violated either be- 
cause of the velocity drift if T 7^ Tq, either because of the 
mismatch in the pulse frequencies. This opens quantum 
channels out of the levitating paths and induce losses in 
the suspended sample. By measuring the sample popu- 
lation after a given number of cycles with resonant fre- 
quencies a;i^2,3,4:='^i234 ^-^d for different values of the 
period T, one can thus determine Tq and measure the 
local gravitational acceleration through Eq.(22). This 



measurement indeed also rests on the determination of 
the recoil velocity Vr — hk/m, but the latter quantity 
is usually known with an excellent accuracy. It is worth 
noticing that the recoil velocity may also be determined 
independently by the present experiment. When the pe- 
riod T is fixed to T := Tg, the suspended cloud looses 
atoms if the pulse frequencies wi,2,3,4 are shifted from 
their resonant values. These mechanism can be exploited 
to lock the lasers performing the suspension on these res- 
onant frequencies. 

We recall here the discussion of Ref.|7j showing that 
the proposed setup is analog to an atomic resonator in 
the momentum picture. It is insightful to consider, as 
plotted on FigjSj the motion of the atomic wave-packets 
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in the energy-momentum picture when both resonance 
conditions are satisfied. The energy includes contribu- 
tions from the internal level energy, from the kinetic en- 
ergy and from the gravitational potential energy. For the 
non-relativistic motion considered here, the energy of an 
atomic packet is a quadratic function of its average mo- 
mentum. The energy-momentum curves are thus a family 
of parabolas, labeled with the internal level and with the 
average altitude of the atomic packet. Each star on the 
Fig|3] stands for an atomic wave-packet. Since the atomic 
motion is conservative and accelerated downwards, these 
stars follow an horizontal leftward trajectory. FigjS] re- 
veals that, in the momentum picture, the atomic motion 
is periodic and occurs in a region bounded by two well- 
defined values ±hk corresponding to the atomic recoil. 
This property strongly suggests an analogy with a mo- 
mentum resonator. Momentum confinement is provided 
here by destructive interferences which shut off the quan- 
tum channels going out of the considered region. 



Total Energy 



^1^^ Total Energy 

(z = z„ + h,l 




FIG. 3: Motion of the atomic wave-packets in the energy- 
momentum picture for the interferometer duration 2To. 
Figs. (a,b,c,d), associated with the 1st, 2nd, 3rd and 4th light 
pulses respectively, show the packets present in coherent su- 
perposition (full stars) immediately after - or transferred 
(transparent stars) during - the considered pulse, whose ef- 
fect is represented by a full red arrow. 



III. SCATTERING ON A SEQUENCE OF 
SHORT RESONANT LIGHT PULSES. 

As a preliminary step toward the prediction of the in- 
terferometer fringe pattern, we compute the macroscopic 
wave-function associated with the levitating atoms in 



a sequence of resonant 7r/2 pulses. The transverse de- 
grees of freedom play a trivial role in the considered 
setup, so that one can follow without loss of general- 
ity a one-dimensional approach. In the following, the 
time-dependence of the momentum width is given by 
Wp{t) = h/y^wl + h?{t — to)'^/wQ. It wiU be simply 
noted Wp. 

We shall use the following property - proven in Ap- 
pendix -. After N pairs of 7r/2-pulses, i.e. at the time 
tN = NT, the wave-function associated with the levitat- 
ing atoms is of the form 



1 

E 



(24) 



Levitating Paths 

The altitudes Zf{^,tpf) are the end points of the tra- 
jectories followed by the atomic wave-packets of maxi- 
mum upward momentum, i.e. which pick up exactly a 
quantum of upward momentum hk during each pair of 
copropagating light pulses. With a fine adjustment of 
the interpulse period T, these packets are maintained in 
suspension, theqrefore we refer to these atomic packets 
as the "levitating packets" and to their trajectories as 
the "levitating paths". We emphasize that the wave- 
function defined by Eq.( |24[ ) is not normalized to unity, 
more precisely its square integral gives the levitating frac- 
tion p\Qy{tN) = / dp|V'lev(^'' tivae In, all 
the levitating wave-packets have acquired the same num- 
ber of upward momentum quanta. However, these pack- 
ets are located at different heights, depending on the 
times at which the successive momentum transfers hap- 
pened. It is convenient to assign to each atomic wave- 
packet a 2iV-component vector ~^ — (ei, e2N) describ- 
ing the history of these momentum transfers. We set 
tj = 1(— 1) if the wave-packet acquired a quantum of up- 
ward (downward) momentum at the j^^i pulse, and = 
if the wave-packet did not go through a transition. The 
levitating atomic packets pick up exactly one quantum of 
upward momentum for each pair of copropagating light 
pulses, i.e. their vector ~^ satisfies 



e2n 



^2n 



1 for 



1, 



.N . 



(25) 



We note Cjs/ the set of such vectors. From the condition 
above, it is obvious that the vectors ~^ associated with 
levitating end points satisfy Vj € {1, 2A^}ej > 0. The end 
point Zf{~^ , tjv) can be expressed as a simple function of 
the vector 



N 

+ ^ [ e2„-iVr {In - {n - l)r) + e2nVr {tN - riT) 



Vr = hk/m is the recoil velocity picked up by the wave- 
packet at each momentum transfer. For the levitating 
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trajectories, one can simplify the expression above by 
using that ~^ G Cj\f, which gives 



with 



Zf{~t,tN)= ^OAT +I]„=i (e2n-l -e2„)wrf ■ (26) 



with 



zm = zo + i^N'^{vr ~ gT)T 



(27) 



In order to compute analytically the sum of Eq.(24), 
we express it as a function the vector ~t . Using Eq.(26l, 
one obtains 



V'lev(?''*Jv) 



p g 



■izoN{p-PN) 



The global phase ^at , depending on the phase convention 
chosen for the lasers, is not relevant to predict the fringe 
pattern of the interferometer. Thus we do not specify 
here its exact expression. We have used the property that 
any vector ~^ £ £_\f satisfies X]^=i ~ The sum in 
the right hand side is analogous to an effective canonical 
partition function in which the configuration would be 
described by the vector Since the coordinates satisfy 
Ej € {0, 1}, and since the pairs of values (e2Tii , e2ni+i) and 
(e2n2 J e2n2+i) ^rc summed independently for ni ^ ri2, 
this sum factorizes. The statistical analogue would be an 
assembly of double wells containing exactly one particle, 
the double wells being independent from one another. 
The two wells have here effective energies that are of 
opposite sign but equal in absolute value. The usual 
dimensionless quantity jSE is replaced by the complex 
number ±Vr{p — pN){iT) /h. Using the fact that each 
double well is independent, one can recast this effective 
partition function as 



1 



■tzojv(p-piv) 



X COS 



N 



VrT 



[P - Pn) 



(29) 



This expression shows that, as announced previously, the 
levitating wave-function experiences an exponential lo- 
calization in momentum space as the number of pulses 
increases. This localization can be interpreted as a fil- 
tering in momentum operated by the multiple-wave in- 
terference. The broadening of the sample size along the 
vertical axis, which happens through the splitting of the 
levitating trajectories, reflects the back-action of such 
momentum localization. The analogy between wave- 
function and partition function is reminiscent of the form 
of the propagator for atomic waves. 

Using Eq.(29), one can evaluate the fraction of the ini- 



tial cloud maintained in levitation after a N pulse pairs. 



= F 



kWn 

— ^T, 2N 
2m 



(30) 



Fi{a,n) 



-|-oo 



due"" cos"(au) (31) 



Clearly F{a, n) ~ 1 only if a ^ 1, which is the condition 
on the coherence length derived earlier. This condition 
must be satisfied in order to obtain interferences between 
two atomic packets flying on different paths. However, 
we shall see that a certain sample fraction can still be 
suspended even if this condition is not met. Levitation 
happens then through constructive interferences at the 
intersection of the different paths, i.e. in the nodes of 
the suspended network. 



IV. MULTIPLE- WAVE ATOMIC LEVITATION 
IN THE SHORT PULSE REGIME. 

In this Section, we compute the fringe pattern of the 
multiple-wave atom interferometer with a simplified 
model for the light pulses. Such model is valid when 
their duration r is much shorter than the resonant 
period Tp, and when the pulses are nearly- resonant, 
i.e. detuned from a frequency Suj ^ l/r, which is 
the regime considered for this setup. The light pulses, 
acting as temporal atomic beam splitters, can then be 
modeled through an instantaneous Rabi matrix with 
elements of equal norm. This model yields excellent 
agreement with the numerical simulations performed in 
the original proposal [7]. In the considered regime, the 
different levitating paths are equally populated, and the 
variation in the total levitating atomic population arises 
purely from interferences effects. In this sense, the short 
pulse limit is the paradigm of the levitation through 
multiple-wave atomic interference and corresponds to a 
perfect quantum trampoline. 



A. Clock fringes. 

Here we still assume that the interpulse-period T is 
fixed at its resonant value T := Tq, but the pulse fre- 
quencies u}j are now shifted from their resonant values. 
As before, the pulse phases are synchronized in order to 
yield a relative phase (j)c = for each pair [See Eq.( 10 1 of 
Sec|l]. In order to understand how these frequency shifts 
affect the levitating wave-function, we consider again 
the scattering of a Gaussian lower-state wave-packet [See 
Eq.(42)], centered around the position Zi, on two upward- 



travelling light pulses interspaced by a duration Tq and of 
respective frequencies uji — uii + 6u)i and LU2 = ^2 ~^ ■ 
We know from Sec|T]that if the frequencies wi^2 were res- 
onant, the path a and b would yield an equal phase pro- 
portional to the invariant / of Eq.( [2l|). Considering the 
expression for the laser phases i n Eq.(|10[ ) and the defini- 
tion of the quantities la.b in Eqs.( 19|20 1, one sees that the 
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frequency shifts Sloi and (5w2 induce a phase imbalance 
between the paths a and b: 



Ia{zi,pN,Ta) = I{k,pN,To) - HSuj2 



Tn 



Tn 



(32) 



Fohowing the hues of Section |Tj the levitating wave- 
function after the pair of upward-traveling pulses is given 
by Eq.(43) with the values above for la^b- Using the re- 
lations Zfa = Zi — hkTo/{2m) and Zfb — Zi -t- hkTo/{2m), 
one obtains 



V'lev&(P'^^'+l) 



„i(*jv+/(fe,pjv,To)/;i) 
2Ar+l^l/4 



X e 

X 



(p-pjv+i) 

^ p-*Tr (p-pn) 



(33) 



\ ^ g4'i( '5'^i-fe(p-PJV + l)/m)-e2(i5'^2- /£(p-pjv+l)/m)] ^ 



We sum only on the values (£1,62) G {0,1}^ such that 
€1 + 62 = 1- One must also consider the scattering of 
a Gaussian upper-state wave-packet, centered initially 
around on a pair of downward-traveling pulses sep- 
arated by the duration Tq, and whose frequencies are 
noted CLI3 and 0-14. Since the gain of an upward momen- 
tum is then associated with a photon emission instead 
of a photon absorption, one simply needs to replace in 
the previous expression the frequencies uji and i02 by the 
values and —1^4 respectively. The result reads 



^leVab'*W+2) 



2^+27ri/4^/^ 



(p-PN+2) 



g-«Tf (p-Pn) 



(34) 
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Similarly we sum on (£3, £4) e {0, 1}^ satisfying £3 -I- £4 = 
1. /~(fc,p7v+i,To) denotes the invariant associated with 
the sample propagation in a pair of resonant downward- 
travelling 7r/2-pulses. Its expression, not needed here, 
can be obtained by replacing the resonant frequencies 
and respectively with the values — cjg = — wSaud 



4 ~ ~uj^ in the expressions of la and lb [Eqs.( T9|20 



Using the Eqs.( 33|34 ), one can show recursively that t 



lie 



wave-function after N = 2M pulse pairs is of the form 



V'leva(P'^w) 



1 



2^7ri/4^/^ 



p (35) 



p 



X C'(-i^,7j-)e-*^/(^'*")(P-P") 



with 



2JV 



As in the previous Section, the factorization of this ex- 
pression is easily obtained, either by using the former 
analogy of a two-particle partition function, either by 
using directly the recursive relations of Eqs.( 33|34 |: 



V'lev(P'^Jv) 



1/4. 



'W. 



p 



/ Suj4n+l + 5uj4n+2 _ Hp - Pn) \ Tq 

\ 2 



m 



M-l 

X cos 

Y SuJ4n+3 + SoJin+A k{p - P n)\ Tq 

V 2 m 



X cos 



(36) 



Again, we do not specify the global phase factor (j)'^. 
We can now describe the clock fringes assuming that all 
pulses are equally detuned, i.e. V i Suji — Suj. The levi- 
tating population yields simply 



Plev 



(37) 



with 
F2(a,n, ( 



due ^ cos" {au -h (f)) cos" {au 



(38) 

We consider here the limit of long sample coherence 
length, which corresponds to a ^ 1, and which is also the 
regime discussed in Ref. 7^. In this limit, the levitating 
fraction becomes simply 



2N I s: To 

Plev = cos ( dwy 



(39) 



This gives readily the full width at half maximum 
(FWHM) Auj, yielding by definition pigy(±iAcj) = 1/2, 
i.e. Aw = 4 arccos (1/2^/(2^)) /Tq. In the limit where 

iV > 1, one easily sees that Aw Tq/Vn. A Tay- 
lor expansion in the equation above yields pYQy{6io) ~ 



Aw 



y21og(2) 



TqVN 



(40) 



One retrieves the usual improvement of the frequency 
resolution Aw with the interrogation time Tjj^^ = NTq 

— 1/2 

as Aw cx , which is typical of atomic clocks. These 

expressions for the width of resonance are plotted in the 
Fig.Q. This figure shows that the width of resonance 
obtained in Eq. ( 40 ) with a Taylor expansion is accurate 



even for a small number of pulses. Excellent agreement is 
obtained between the expression of Eq. ( 39 ) and the band- 



width computed from the simulations used in Ref. . As 
in Ref. [3, we considered the mass m of ^^Rb atoms and 
a wave- vector k ~ 10 ''m"^. 
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FIG. 4: Full width at half maximum (KHz) for the levitating 
population as a function of the number of pulse pairs. Full 
red curve shows the bandwidth obtained with Eq.([39|. Blue 



dashed curve denotes the simplified expression of Eq.(40l 



The dots denote the bandwidth obtained from the numeri- 
cal simulations of the levitating wave-packets used in Ref.[7]. 



B. Gravimeter fringes. 

When the setup operates as a gravimeter, one scans the 
interpulse period T around its resonant value Tq. The fre- 
quencies are chosen as resonant for an average atomic mo- 
mentum p'q at the beginning of the four-pulse sequence, 
i.e. uj4j+i = LOij+i = and uJij+2 = ^4^+3 = UJ2, with 



UJ12 satisfying Eq.(16) with p'q. We note pj is the mo- 
mentum of the levitating wave-packets at the the time 
tj = jT. The mismatch in the resonance condition in- 
duces a drift in momentum, i.e. pj = po + j{hk — mgT). 
One sees on the resonance condition [Eq.(16)] that the 



mismatch of the atomic momentum with respect to the 
resonant value p'^ is responsible for the frequency detun- 
ing of the jth pulse 

= (-1)^"*"^ - k - Po + C(j)(^fc - mgT) 
ml 

The first factor in the right hand side accounts for the re- 
versal of the direction of the wave- vector k between each 
pair of pulses. We have introduced for convenience the 
auxiliary function C,{j)T which gives the instant at which 
the j^Yi pulse is performed in the considered sequence: 
C(4m-hl) = 2m, C(4m-|-2) = 2m -hi, C(4m-h3) = 2m+l 
and C(4m -f 4) 2m + 2. The previous discussion on the 
path imbalance induced by a detuning can be readily 
applied. One obtains a levitating wave-function given by 
Eq.(41 1 with the frequency detunings 5ujj replaced by the 
values of the equation above. In the limit of long con- 
densate coherence length, the levitating fraction reads 



Plev 



2M 

[]C0S2 

i=i 



Losses in the levitating cloud arise from the initial mo- 
mentum shift and from the mismatch in the resonant 



period. The gravitational acceleration g is determined 
from Eq.(22). One assumes that the recoil velocity is 



known with an excellent precision, so that the precision 
of the measurement of the acceleration g is simply lim- 
ited by the determination of Tq, i.e. /S.g/g — ATq/Tq. 
Considering that a fractional loss e can be detected in 
the levitating atomic population, one can determine the 
resonant period Tq up to an uncertainty ATq = ^2 — Ti , 
with Ti and T2 such that Pa{Ti) = Pa{T2) = 1 — e. Tak- 
ing the logarithm of Eq.( |4ip , using a Taylor expansion 
and replacing the sum by an integral, one finds 
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v/-61og(l-e) 1 



kgT§ 



(41) 



As in Ref. [51 [H] j the atomic velocity drift over the succes- 
sive cycles yields a measurement accuracy scaling with 
the number N of cycles as Ag/g c>c l/N^^"^. The ac- 
curacy given above agrees very well with the accuracy 
predicted by numerical simulations of the sample. Both 
are plotted on the following figure. 

Ag/g 
10 



10-2 
10-3 
10-4 
10-5 



10 

N 



100 



FIG. 5: Relative uncertainty on the gravitational acceleration 
predicted by Eq.(41| (red dashed curve) and that obtained 
with numerical simulations (blue curve) as a function of the 
number of cycles A'' (corresponding to 2N pulse pairs). 



V. EXPERIMENTAL REALIZATION OF THE 
SETUP. 

We investigate here suitable atomic species for the con- 
sidered scheme. The considered two-level atom should 
have a clock transition whose life-time should fulfill two 
seemingly contradictory requirements. In order to in- 
terrogate the atoms over several bounces, this life-time 
should be much longer than the interpulse period Tq. 
It should also be short enough as to enable the realiza- 
tion of a 7r/2-light pulse in a duration less than Tq/2 
with a realistic laser intensity. These two criteria imply 
a non-trivial trade-off between sharpness linewidth and 
coupling on the transition. 

Besides, the atomic sample should be cold enough so 
that the short 7r/2 light pulses address its entire velocity 
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width. In order to preserve the full atomic population in 
levitation, one should use an atomic sample with a coher- 
ence length much longer than the separation between two 
paths. As discussed here-after, this requirement is not 
met by fermionic samples: a Bose-Einstein condensate is 
required to obtain such global interferences. Neverthe- 
less, partial levitation can still be achieved with fermions 
for a long duration. 



A. Local interferences with Ultra-cold Fermions. 

The compromise on the clock transition is ideally 
achieved by the ^Sq — '^Pq transition of the fermionic iso- 
tope ^^-"^Yb TS], which has a wave-length A = 578.4 nm 
and a frequency width of 10 mHz. For this transition, a 
7r/2 pulse of 4/is (on the order of Tq/IOO) can be realized 
with a laser power of 0.3 W focalized on a spot of radius 
1 mm (1/e^ beam radius) |15j . 

Even if temperatures as low as O = 130 nK have been 
achieved experimentally by sympathetic cooling with a 
cloud of 10^ atoms of this isotope [TB], well within the 
regime of quantum degeneracy, the associated dispersion 
in momentum is unfortunately too wide to obtain inter- 
ferences between atomic packets propagating on differ- 
ent paths. Indeed, the coherence length w is still much 
smaller than the distance L between two adjacent paths, 
i.e. L/w — hk^y/ksQ/m'^^^ ^ 10 with ks the Boltzmann 
constant. Nevertheless, the interferences of the wave- 
packets at the nodes of the network still enable to levitate 
a significant fraction of the atomic population and yield 
a sufficient contrast. Indeed, one can tailor either con- 
structive, either destructive interferences at these nodes. 
We note that the interferences at each node correspond to 
those of an atomic Borde-Ramsey interferometer issued 
from a node located at the same height and occurring at a 
time 2T earlier. When {Ib — Ia)/fi = 0(7r), the arms of the 
Borde-Ramsey interferometer yield the same phase and 
thus constructive interferences at each node. In contrast, 
when {Ii, — Ia)/h = 7r/2(7r) , one obtains a phase detuning 
of TT between the arms of the Ramsey-Borde interferome- 
ters yielding destructive interferences at each node where 
several trajectories intersect. These interferences quickly 
shut off the levitation as the number of pulse pairs in- 
creases: beyond ten pairs of pulse, the levitating frac- 
tion is found below 10~^, so that no atom in the sample 
is expected to levitate. The contrast is defined as C = 
(Plev(0)-?'lev('^/2))/(Plev(0)+Plev('^/2)): where pigy(O) 
and pYQy{'K/2) are the levitating fraction for constructive 
and destructive interferences respectively. Fig j6] gives the 
levitating fraction at resonance as well as the contrast of 
the interference pattern for a sample of ultra-cold ^^^Yb 
atoms. The resonant period is niY\, — 0.4 ms for this 
isotope. Starting with a sample of lO'' atoms, one ex- 
pects to maintain roughly 400 atoms in suspension after 
100 pairs of pulses which is still enough to be detected 
with state-of-the-art techniques. 




°1 10 20 30 40 50 60 70 80 90 100 



FIG. 6: Levitating fraction (red curve, left scale) and contrast 
(blue curve, right scale) as a function of the number of pulse 
pairs for a sample of Yh atoms at a temperature = 130 nK. 



B. Global interferences with ultra-cold Bosons. 

In order to obtain interferences between atomic pack- 
ets flying on different paths, the initial sample must have 
a sufficient coherence length. Such coherence is only 
achieved by BECs, so that this time bosonic species need 
to be considered. In order for the setup to operate as a 
clock, the considered atom should own a narrow transi- 
tion in its internal structure. For this purpose, one can 
consider using the forbidden optical transition ^Sq— '^Po, 
enhanced by a dc magnetic field, of neutral even isotopes 
of Yb, Sr, Mg, or Ca [T7]. In the following discussion, we 
use the atomic transition parameters presented in this 
reference. 

We focus on the two species ^^"^Yb and '^''Ca, which 
have been Bose-condensed recently [TB], and which have 
different advantages. Since both the recoil velocity and 
the resonant period Tq decrease with the atomic mass, 
the requirement on the coherence length is less stringent 
for heavy atoms. This point advocates for the use of 
heavy ultra-cold samples such as the Yb isotopes, for 
which the coherence length required to obtain inter-path 
interferences is roughly — 1 A*™., to be compared 

with vjQg^ ~ 10 /xm for Ca atoms. On the other hand, 
since ^ "^o Yb' ^ setup using ^°Ca atoms involves 

a smaller Rabi pulsation on the clock transition. This 
property benefits to Ca atoms, since a smaller Rabi pul- 
sation requires less intense magnetic and light field to per- 
form the transition, and thus smaller spurious frequency 
shifts associated with the second-order Zeeman effect and 
the Stark shift respectively. We impose a pulse duration 
T = 0.3To for the 7r/2-light pulses, which corresponds to 
an operation of the multiple-wave interferometer in the 
long-pulse regime, more precisely with 60% of illumina- 
tion time. We assume that the dc magnetic field B can 
be controlled up to a relative accuracy of 10^^, while the 
light intensity / is monitored with a relative accuracy 
of 10"**. Using ^^''Yb atoms, one performs the desired 
light pulses by applying a magnetic field B — 200 G and 
a laser power P ~ 50 mW focused as above on a spot 
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of radius 0.1 mm, yielding a magnetic frequency uncer- 
tainty A^Yb ~ 0.50 Hz and a light shift uncertainty 



0.46 Hz. With '*°Ca atoms, one realizes instead 
the pulse with a smaller magnetic field i? = 40 G and a 
laser power P ~ 120mW focused on a similar spot, yield- 
ing the respective magnetic and light shift uncertainties 
Ag ~ ~ 0.27 Hz. These figures can seem high 

in comparison with other schemes performing an interro- 
gation with a Rabi pulsation in the Hz range such as op- 
tical network setups. These significant shifts and uncer- 
tainties are caused by the short pulse constraint imposed 
by the value of the resonant period Tq. This constraint 
would be largely relaxed when operating this architec- 
ture with an atomic transition in the UV range |19| : the 
period Tq = (h/mc) x {oj/ g) would become much longer, 
thereby enabling atomic levitation with weak laser fields 
inducing smaller frequency shifts. 



VI. CONCLUSION. 

It is appropriate to compare briefly the theoretical 
method used here with the treatment of other recent 
atomic levitation experiments [5", '11] . These experiments 
rely on the illumination with standing light waves instead 
of the traveling waves considered here. 

In all the atomic levitation experiments [51[7|[SllU'j the 
analysis of the sample sustentation relies on the articula- 
tion between two distinct regimes of propagation, namely 
the atomic propagation in the dark and the atomic prop- 
agation in the laser field. The propagation in the dark 
has been solved exactly with the atom optical ABCD for- 
malism |12j under the following assumptions: the atoms 
evolve in an external potential which is quadratic, and 
atomic interactions are sufficiently weak to be neglected. 
Such assumptions are legitimate for typical clock exper- 
iments. In the absence of illumination, atoms are sub- 
ject only to the local gravitational potential which is 
very well-approximated by a quadratic potential. Mean- 
field effects can be included in the ABCD method for 
moderate interactions when a single atomic cloud is in- 
volved [50], or if the system is loaded with an atom 
laser [3T] . So far our treatment is equivalent to that used 
in the Refs.[9l[IT]. 

However, a clear distinction needs to be made between 
our proposal and the experiments [SJ [IT| regarding the 
propagation in the laser field. During the sample illu- 
mination, the light field couples plane atomic waves of 
different momenta. In our proposals involving travel- 
ing waves, the light field couples each atomic momentum 
eigenstate |p) to a single momentum eigenstate |p -I- /ik) 
as expected from the atomic recoil. Such scattering, in- 
volving a single quantum channel, can be effectively mod- 
eled by means of a 2 x 2 Rabi matrix. In contrast, the 
scattering resulting from the illumination with a stand- 
ing wave involve multiple quantum channels and couple 
instead each atomic momentum eigenstate |p) to a mul- 
tiplicity of states {|p-t-m/ik), m € Z}. The description of 



the scattering is thus much more complex when standing 
light waves are involved instead of traveling waves. It is 
not clear whether a similar summation with an effective 
partition function can be obtained when standing waves 
are used. One should note, however, that in practice for 
the experiments reported in the Refs.[9l [11], the pulses 
are tailored so that the standing wave scattering involve 
predominantly the first few orders - in Ref.[ll|. scatter- 
ing beyond the second order can be safely neglected -. 
A finite-dimensional Rabi matrix is then also suitable do 
describe the atomic scattering by standing light waves. 

To summarize, we have developed a theory of atom 
interferometers using multiple-wave levitation. We have 
used a propagation invariant to describe the atomic scat- 
tering by a periodic series of short light pulses. The sensi- 
tivity of the sensor as an atomic clock or as a gravimeter 
has been computed, and yields excellent agreement with 
the numerical simulations presented in Ref. [7]. We have 
proposed to use either fermionic or bosonic species to 
explore localized or global interferences in the levitating 
atomic sample. Using recent advances in the cooling and 
spectroscopy of Yb and Ca isotopes, we have shown that 
such clock system can be implemented with state-of-the 
art experimental techniques. Further developments will 
include a more sophisticated treatment of atomic beam- 
splitters in order to explore the system behavior in the 
long pulse regime. Last, we mention that our scheme 
would be particularily well-suited for clocks running with 
high-frequency transitions in the UV range |19| . 
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Appendix: Expression of the levitating wave-packet. 



We establish by recursion the expression of Eq. ( 24 ) 

(P"PJV)^ 



1 



E 



= <^ •tN)ip-PN) 



Levitating Paths 

For = 0, the sum in this equation contains a sin- 
gle term which is the initial lower-state Gaussian wave- 
packet expressed in the momentum picture. This proves 
the desired property at the rank A^ = 0. Let ut as- 
sume that this property is satisfied at the rank N. To 
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establish it at the rank N + 1, we consider the propa- 



gation of the wave-function of Eq.(24) in a sequence of 
two hght pulses of short duration, of respective initial 
instants ijv — NT and ^at+i — {N + 1)T. We assume 
that the levitating wave-packet is in the lower state at 
the instant t — t^, and that the light pulses are upward 
travelling. Since the propagation is linear, one can con- 
sider independently each Gaussian wave-packet present 
in the sum of Eq.(24). It is thus sufficient to prove the 



two following facts. Any Gaussian packet in the sum of 
Eq.(24|, associated with the altitude z^, is split into two 



levitating wave-packets, whose central altitudes Zfa and 
Zfh are obtained by following either the lower or the up- 
per arche issued from the altitude Zi [See Fig{l]. Besides, 
these two wave-packets, written in the form of Eq.p4|), 
bear a phase (fiN+i which is independent of the altitude 
Zi of the initial wave-packet. 

Following these lines, we consider the following wave- 
function at the instant tjq 







^ 't2^i(p-Piv) / 

(42) 

The lower-state wave-function is a Gaussian wave- 
packet in the momentum picture, which would be noted 
wp{p,Zi,pN,w,tN) with the conventions of Eq.([2|. To 
compute its propagation between the instants In and 
tAr_|_i, one can apply the discussion of Section |lj which 
gives readily the levitating wave-packet at time tjv+i- We 
discard here the lower-state component of Eq.(|14|), since 



the corresponding wave-packets have missed one quan- 
tum of upward momentum and thus do not contribute 
to the levitating wave-function. Using Eqs.( |13|r8 1, this 
wave-function can be expressed with the path invariants 
noted la = Ia{zi,pN,T) and h = Ib[z.i,pN,T) 



J<S>N 



(p-PjV + l) 



V'lev.(P'^^+i)=2A^+i^i/4^e 



(43) 



Pat -1-1 = Pn + hk ~ mgT is the momentum of the levi- 
tating packets at time t = tAr+i. As outlined in Section 
m when the pulse frequencies satisfy the resonance 
conditions, the quantities and /(, are equal and given 
by Eq.(21|, i.e. la = h = I(k,PN ,T). One then obtains 
a phase factor (/)7v+i — 4>n + I{k,PN, T) -f 7r/2 which is 
independent of the altitude z^, and consequently equal 
for all the wave-packets present in the sum of Eq.(24|. 



Besides, when the altitude Zi takes all possible values of 
the end points z/ of the levitating trajectories at time 
t = tpf, one obtains with Zfa and Zfb all the possible 
end points at the time t — ^at+i. The desired property 
is thus proven at the rank N + I. Should we have con- 
sidered instead the transfer of atomic wave-packets from 
the upper to the lower state with downward-travelling 
pulses, the argument would have been naturally identical. 
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